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NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 

TECHNICAL  NOTE  2159 

ON  THE  PARTICULAR  INTEGRALS  OF  THE  PRANDTL  -BUSEMANN 
ITERATION  EQUATIONS  FOR  THE  FLOW 
OF  A  COMPRESSIBLE  FLUID 
By  Carl  Kaplan 

SUMMARY 


The  particular  integrals  of  the  second-order  and  third-order  Prandtl- 
Busemann  iteration  equations  for  the  flow  of  a  compressible  fluid  are 
obtained  by  means  of  the  method  in  which  the  complex  conjugate  variables 
and  z  are  utilized  as  the  independent  variables  of  the  analysis.  The 
assumption  is  made  that  the  Prandtl-Glauert  solution  of  the  linearized  or 
first-order  iteration  equation  for  the  two-dimensional  flow  of  a  com¬ 
pressible  fluid  is  known.  The  forms  of  the  particular  integrals,  derived 
for  subsonic  flow,  are  readily  adapted  to  supersonic  flows  with  only  a 
change  in  sign  of  one  of  the  parameters  of  the  problem. 


INTRODUCTION 


For  the  past  several  years  iteration  methods  have  been  increasingly 
applied  to  the  solution  of  compressible-flow  problems.  The  most  useful 
one  from  the  point  of  view  of  aeronautical  applications  and  the  one 
discussed  in  this  paper  is  based  on  small  perturbations  with  respect  to 
the  undisturbed  flow.  The  Prandtl-Glauert  and  Ackeret  solutions  in  two- 
dimensional  subsonic  and  supersonic  flow,  respectively,  obtained  by 
means  of  the  linearization  of  the  fundamental  nonlinear  differential 
equation  for  compressible  flow,  are  presumed  to  be  known  and  are  taken 
as  the  initial  steps  in  this  iteration  process.  Higher-order  solutions 
are  then  obtained  by  retaining  appropriate  powers  and  products  of  the 
perturbation  quantities.  This  method  of  iteration  has  been  variously 
labeled  the  Ackeret  iteration  process  and  the  Prandtl-Busemann  small 
perturbation  method  when  limited  to  two-dimensional  subsonic  flow.  The 
procedure  has  been  extended  in  recent  years  to  both  two-dimensional  and 
axisymmetrical  supersonic-flow  problems. 

In  a  recent  publication  (reference  l).  Van  Dyke  succeeded  in 
obtaining  by  trial  the  particular  integral  of  the  nonhomogeneous  second- 
order  iteration  equation  for  the  velocity  potential  in  supersonic  flow. 
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The  general  solution  is  then  easily  obtained  by  adding  solutions  of  the 
homogeneous  equation  with  proper  regard  to  the  boundary  conditions  at 
the  surface  of  the  solid  and  at  infinity. 

The  purpose  of  the  present  paper  is  to  show  a  procedure  by  means 
of  which  the  particular  integrals  of  the  higher-than-f irst- order  iteration 
equations  can  be  derived  in  a  systematic  manner.  The  explicit  expression 
obtained  for  the  particular  integrals  of  the  second- and  third-order 
iteration  equations  are  believed  to  yield  essentially  the  solution  of 
the  problem  of  high  subsonic  flow  past  an  arbitrary  two-dimensional 
profile,  since  it  is  never  a  difficult  problem  to  supply  the  solutions 
of  the  homogeneous  equation  necessary  for  the  fulfillment  of  the  boundary 
conditions.  It  is  noteworthy  that  the  particular  integrals,  derived  for 
subsonic  flow,  can  be  adapted  to  supersonic  flow  with  simply  a  change 
in  sign  of  one  of  the  parameters. 

FUNDAMENTAL  EQUATIONS 


The  fundamental  nonlinear  differential  equation  governing  the  flow 
of  a  compressible  fluid  is 

(c2  .  u2)|  +  (c2  .  v3)|  -  uv(|  +  |)  -  o  (1) 

where 

X,  Y  rectangular  Cartesian  coordinates  in  flow  plane 

u,  v  fluid  velocity  comppnents  along  X-  and  Y-axis,  respectively 

c  local  speed  of  sound 

The  condition  for  irrotational  motion  is  that 


Su  _  Sv 

By  "  Bx 

and  leads  to  a  velocity  potential  $  defined  by 
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If  the  body  is  held  fixed  in  a  uniform  stream  of  velocity  U,  the 
relation  between  the  local  speed  of  sound  c  and  the  speed  of  the 

fluid  \Ju^  +  is  given  for  adiabatic  processes  by 


6  -  1  +  Wf  - 

where 

sound  speed  In  undisturbed  fluid 

7  ratio  of  specific  heats  at  constant  pressure  and  constant  volume 
Mach  number  of  undisturbed  stream  (U/cot) 

With  the  introduction  of  a  characteristic  length  Z  as  unit  of 
length  and  the  undisturbed  stream  velocity  U  as  unit  of  velocity,  the 
quantities  X,  Y,  u,  v,  and  $  for  the  remainder  of  the  analysis 
denote,  respectively,  the  nondimens ional  quantities  X/Z,  Y/Z,  u/U, 
v/U,  and  $/UZ,  while  c  and  c^  retain  their  original  meanings.  By 
means  of  equation  (2),  equations  (l)  and  (3)  then  become,  respectively. 


and 


2mJVy*XY  =  ^ 


(*0 


1  + 


2L - 
2 


hi  2 


(u2  +  v2) 


(5) 


where  the  subscripts  X  and  Y  denote  partial  differentiations  with 
respect  to  the  designated  variables. 


In  order  to  obtain  the  iteration  equations  based  on  small  perturbations 
of  the  undisturbed  stream,  the  assumption  is  made  that  the  velocity 
potential  $  can  be  expanded  in  the  form 


0  =  X  +  ^  +  $2  4-  $2  +  . . . 


(6) 


4. 
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For  the  purpose  of  defining  and  controlling  the  iteration  procedure, 
the  function  $n+1  and  its  derivatives  are  then  regarded  as  small 

compared  with  the  preceding  approximation  $G  and  its  derivatives. 
From  equations  (2)  and  (6), 

u  =  =  1  +  +  °2X  +  ^3X  +  •  •  • 


and 


Y  =  Gy  =  +  $2Y  +  ^3Y  +  •  •  • 

When  these  expressions  for  u,  v  are  introduced  into  equation  (4), 
together  with  the  expression  for  c^/c^  given  by  equation  (5),  and 
the  powers  and  products  of  ^  and  their  derivatives  are  grouped 
according  to  the  assumptions  of  the  small  perturbation  method,  the 
following  iteration  equations  for  the  first  three  approximations  $]_, 
and  $3  result: 


i1  ■  M.2)*1XX  *  *LYY  ■  0 


(7) 


1  -  H»2Wxx  +  ^YY  =  M2  (7  +  1)°1XC’1XX  +  (7  -  1)^>ix^iYY  +  ^lY^lXY 


(8) 


^1  -  ®WV  +  =  M  + 


3XX  ‘r  '"3YY  =  +  °2x)  <7  *  U^IXX  +  (7  ' 


^lY2^  -  iHlXX  +  (7  +  1)®1yyJ  + .  ®1X  7  +  1 ) ^2XX 

(7  -  1 )  °2Yy]  +  2(hx$lY$lXY  +  ^2Y®1XY  +  ^IY^Xy) 

✓ 

(9) 


For  slender  bodies,  the  first  few  steps  of  this  iteration  process  may  be 
expected  to  yield  an  accurate  result  with  the  exception  of  a  small  region 
in  the  neighborhood  of  a  stagnation  point.  Even  at  stagnation  points 
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the  iteration  method  has  been  shown  to  represent  correctly  the  effect  of 
compressibility  (reference  2).  The  accuracy  of  the  calculations 
obviously  depends  upon  the  number  of  terms  determined,  each  additional 
term  reducing  the  region  of  inaccuracy  in  the  neighborhood  of  a  stagnation 
point. 

The  iteration  equations  (7)>  (8),  and  (9)  may  be  put  into  more 
familiar  forms  by  the  introduction  of  a  new  set  of  independent  variables  x 
and  y,  where 


Thus,  for  Moo  <  1,  equation  (7)  is  transformed  into  a  Laplace  equation, 
whereas  equations  (8)  and  (9)  are  transformed  into  Poisson  equations 
with  the  right-hand  sides  composed  of,  respectively,  double  products  and 
triple  products  of  previously  determined  perturbation  quantities.  It 
is  further  assumed  that  the  solution  of  equation  (7)  is  available.  This 
initial  step  in  the  approximation  to  the  exact  nonlinear  solution  is 
usually  easily  obtained,  as  it  represents  the  Prandtl-Glauert  approximation 
(reference  3,  appendix  B).  The  purpose  of  the  present  paper  is  then  to 
derive  explicit  expressions  for  the  particular  solutions  of  the  second - 
and  third-order  iteration  equations  (8)  and  (9). 


CALCULATION  OF  THE  PARTICULAR  INTEGRAL  OF  THE  SECOND-ORDER 

ITERATION  EQUATION 


By  introducing  the  independent  variables  x  and  y  defined  by 
equation  (10),  the  second-order  iteration  equation  (8)  becomes 

$2xx  +  $2 yy  =  ^cc?  +  a)^Lx^lxx  +  ^ly^lxy  (H) 


cr  =  Z  +  1  M°° 

2  p2 

P2  =  1  -  M^2 


where 
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and  where  use  has  been  made  of  Laplace's  equation 


°ixx  +  ^lyy  "  0 


The  procedure  for  obtaining  the  particular  integrals  of  the  higher- 
order  iteration  equations  is  based  on  the  use  of  the  complex  conjugate 
variables  z  and  z  as  independent  variables.  Thus, 


z  =  x  +  iy 

z  =  x  -  iy 

and  the  equivalence  of  operators 

a  &  a 

Sc  c)z  +  M 

a  _  Jb  _  a  \ 

^  {&.  5zJ 

+£_ 

OX2  OZ2  5531  OZ2 

£_  _  bi_  ,  a2 

ay2  lzs  *  -  a-2 

Then  equation  (7)  for  0-j.  becomes 

=  ° 

The  most  general  real  solution  of  this  equation  is 


«,  =  i 
1  2 


w1(z) 


wi(z) 


♦ 


* 


(12) 


(13) 


W 
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=  R.P.w-^z)  =  R.P.w^z) 

where  wq(z)  is  an  arbitrary  analytic  function  of  z,  wq(z)  is  its 
conjugate  complex,  and  where  the  symbol  R.P.  stands  for  "real  part  of"* 
The  imaginary  part  of  wq(z)  is  a  function  say,  related  to 

by  means  of  the  Cauchy -Riemann  equations  and  hence  also  satisfies  Laplace1  s 
equation.  The  function  i|/]_  does  not  represent  the  stream  function  of 

the  actual  compressible  flow  and  does  not  appear  in  the  final  expressions 
of  the  particular  integrals.  The  following  relations  will  be  found  use¬ 
ful  and  are  easily  verified: 


*1 

=  |(wx  + 

wx) 

* 

CSI 

I — 1 

hTcvi 

II 

vz) 

* 

4iy 

=  |(vlz  ■ 

-  Wlz) 

^lxy 

2^  vlzz 

-  niz) 

$lxx 

~  2^W1ZZ 

+  ^lzz) 

Then 

^x'flxx  =  R.P.  (viz  +  wiz)wlzz 

2<*ly<Ilxy  =  “  R*P.(wlz  -  ^lz)wlzz 
and  equation  (ll)  for  the  second  approximation  $2  becomes 

°2zi  =  |mm2  R.P.[ctv1zw1zz  +  (2  +  a)wlzwlzzj 

=  K2  R.P.g(wlz2)z  +  (2  +  cr)(wlzwlz)zJ 
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If  $2  is  defined  to  be  the  real  part  of  a  nonanalytic  function  vq  {z,z), 
then 


+  (2  +  a) (v±2?±zj  z 


(14) 


This  equation  can  be  integrated  immediately  by  inspection  and  yields 
the  general  solution 


\ 


w2  =  Jh>o2 


a-  2  , 
2zwlz  + 


(2  +  ffJVlz 


(15) 


where,  because  only  the  real  part  $2  w2  interest,  only  one 

arbitrary  analytic  function  F(z)  need  be  included.  The  function  F(z) 
satisfies  Laplace* s  equation  and  is  so  chosen  as  to  satisfy  the  required 
boundary  conditions  at  the  surface  of  the  body  and  at  infinity.  The 
part  of  the  expression  on  the  right-hand  side  of  equation  (15),  excluding 
the  arbitrary  function  F(z),  is  the  particular  integral  of  equation  (l4) 
and  may  be  expressed  in  real  form  in  the  following  manner: 

Suppo se 


F(z)  =  -  |zwlz2  +  (2  +  a)w1wlz  +  f(z)  (16) 


where  f(z)  is  again  an  arbitrary  analytic  function  of  z. 
Then  with  the  aid  of  the  relation 

”lz  '  §(lr  ■  ‘I?)**1  *  M,l)  ■  '  “ly 

where  use  has  been  made  of  the  Cauchy-Riemann  conditions 


0 

0 


lx 

iy 


* 


iy 


=  -  ♦ 


lx 
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the  expression  for  <f>2>  obtained  from  equation  (15),  becomes 


$2  =  ^M^2  "  2cry^lx(I'ly  +  2(2  +  ff)®l®lx  +  R-P-f(z) 


(IT) 


The  expression  on  the  right-hand  side  of  this  equation,  excluding 
R.P.  f(z),  namely. 


iy 


lx 


(18) 


corresponds  precisely  to  the  particular  integral  obtained  by  Van  Dyke 
(reference  1)  for  two-dimensional  supersonic  flow  with  a  replaced 
by  -a,  where  for  supersonic  flow  the  definition  of  a  is 


a 


7  +  1 


M. 


1 


jM  It  is  rather  noteworthy  that  the  particular  integral  of  the  second-order 

iteration  equation  (8)  can  be  obtained  for  both  subsonic  and  supersonic 
flows  by  simply  interchanging  the  sign  of  the  parameter  a. 


CALCULATION  OF  THE  PARTICULAR  INTEGRAL  OF  THE  THIRD -ORDER 

ITERATION  EQUATION 


In  this  section,  the  particular  integral  of  equation  (9)  involving 
only  and  their  derivatives  is  derived.  For  this  purpose, 

the  variables  x,y  and  the  parameter  a  are  introduced.  Equation  (9) 
then  takes  the  following  form: 


®  3xx  +  ®  3yy  ^oo 

(1  +  a) 


j(l  +  o ) (<i>  1x3!? 2x  +  $lx^2xx)  +  ^2P2(l  +  cr )  -  1 

^lxx®ix2  +  I 


®  lx^!  ly®  lxy  + 


2p2(i  +  a)  -  | 


2  +  ^{crp2  -  1)<1>lxx$ly2  + 


4>lxy®2y  +  <tly3>2xy 


(19) 
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Use  is  again  made  of  the  complex  conjugate  variables  z  and  z  as 
independent  variables.  Thus, 

$2  =  R.P.w2(z,z) 

^2x  =  r*p*(w2z  +  w2z) 

®2y  =  R.P.i(v2z  -  w25) 

^2  xx  =  R*p*(w2zz  +  ^w2zz  +  m2.zt) 

$2xy  =  R.P.i(w2zz  -  Wggg) 

$lxx^lx2  =  R-p-|(wlz  +  ^)2vlzz 
^lxx^ly2  =  -R*P-|(wlZ  -  ^lz)2Wlzz 
®lxx®2x  +  *lx*2xx  =  R-P-!(Wlzz  +  ^1ZZ)  (w2z  +  w2z)  + 

R.p.|(wiz  +  wlz)(w2zz  +  2v2zz  +  w2zz) 

^lxy^2y  =  -R-p-|(wlzz  "  *lzS)(w2z  “  v2z) 

^ly®2xy  =  “  R,p,2(wlz  “  ^lz^w2zz  ”  w2zz^ 
®lx®ly®Ixy  =  -R*P-KW^2  “  ^lz9Wpzz 
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Then  equation  (19)  *  with  §3  =  R.P.w^z^z),  can  be  written  as  follows 


w3zz  =  TfC1  -  P2)<j(w2zw1zz  +  w2zzw1z  +  w2i*Uz  +  w2zzhz)  + 


t-(1  -  P2)(2  +  a)(v_  w, —  +  w0  v.,-  +  w 0-vn  +  wn — w..  ) 
4'  ''  2z  lzz  2zz  lz  2z  lzz  2zz  lz' 


(1  -  P2)a  p2(l  +  2a)  -  (1  -  2a)] 


wlz  wlzz  + 


(1  -  p2)  p2(4  +  5a  +  2a2)  +  a(3  +  2a 


wlz  wlzz  + 


|(1  -  P2)[p2(2  +  5a  +  2a2)  +  (-2  +  a  +  2a2)]  wlzv12vlzz 


Introducing  the  expressions  for  w2(zjz)  and  its  derivatives  with 
respect  to  z  and  z  from  equation  (15)  yields  for  w^g  the  following: 

v3zz  =  T^(1  "  P2)2a2i(vlz2vlzz)z  +  q^1  "  P2)2o(2  +  a)(iw1)2(wlzwlzz)z  + 


3^(1  -  P2)2^(w1zFz)z  +  3^(1  -  P2)aj]p2(2  +  a)  +  (4  +  3a j]  wlz(w 


lz  }z 


..  q2\2 


(1  -  p2rcx(2  +  a)(w1,72)_(zw1 -)-  + 


1^(1  -  P2)2(2  +  0)2wlzz(hhz)g  + 


hi1 


-  P2)2(2  +  <J)(^lzFz)zz  +  -  P2)2(2  +  ff)2(^l2)zwlzzz  + 


-  P2)a2(5  +  3P2)(w1z3)z  + 


-  P2)[p2(8  +  10a  +  3a2)  +■  a(6  +  5a)lv12gwlz2 
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In  the  derivation  of  this  expression  for  w^zg,  free  use  was  made 
of  the  fact  that  insofar  as  the  real  part  $3  of  W3  is  concerned, 
terms  of  the  nature  g(z)E(z)  and  g(z)h(z)  are  equivalent.  It  is 
important  to  note  that  this  type  of  operation  leaves  unaltered  the  real 
part  $2  of  W3.  Since  $3  is  the  quantity  sought  for  in  the  calcu¬ 
lations,  changes  in  the  imaginary  part  of  W3  are  of  no  consequence 
in  the  final  results. 


Equation  (20)  can  be  integrated  immediately  by  inspection  and  yields 
the  following  result: 


1 — 1 

H|co 

11 

-  32)2tf2z2vlz2vlzz 

+  m 

(1  -  p: 

m(1 

-  P2)2ct(2  + 

a)™lVlZVlZZ  + 

1 — 1 

*Ts 

-  p2)a[p2(2 

+  a)  + 

(b  + 

30  )J  V; 

1 — l 

rH|CVi 

-  P2)2a(2  + 

°)Klzwlz2  • 

+  m(1 

TZ^1 

-  P2)2(2  +  c 

>)(«!£ 

+  F  ' 
z 

&P  +  ■ 

^(1 

-  P2)ct2(3P2 

+  5)zw1z3  + 

32^ 

c(6  + 

5°)]  *15  f  wlz2  dz 

2-  2 


(21) 


Equation  (21)  is  the  particular  integral  of  equation  (20).  The 
most  general  solution  is  obtained  by  adding  an  arbitrary  analytic 
function  G(z),  satisfying  the  homogeneous  or  Laplace  equation  Gzz  -  0. 

An  arbitrary  function  of  z,  customarily  included  in  the  general  solution, 
need  not  be  considered  here  because  only  the  real  part  $ ^  of  W3  is 

of  interest^  In  fact,  the  omitted  arbitrary  function  is  the  complex 
conjugate  G(z). 
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In  order  to  obtain  the  desired  form  of  $3  (  tiie  real  part  of  W3) 

from  equation  (21),  F(z)  is  replaced  by  its  expression  given  in 
equation  (l 6),  and  the  real  part  of  f(z)  in  equation  (lj)  is  replaced 

by 


R.P.f(z)  -  +  2ay®lx®ly  “  2(2  + 

M 

00 

The  final  form  of  the  particular  integral  of  the  third-order  iteration 
equation  (9)  then  "becomes 


$3  =  -  |(l  -  32)cry(^2y  +  $ly*2x)  +  |(l  -  P2)(2  +  a)(^2  +  9^)  + 

B(1  ~  -  <tly2)  -  Jd  -  P2)2ff2y2®ix$ly^y  + 


^(1  -  P2)a[(6  +  5a)  +  3(-2  +  a)p2]y<i>ly3  - 
55(1  -  p2)a[(l0  -  a)  -  (10  +  7a)P2]y%x2$: 


^(1  -  P2)2a(2  +  ajy*^ 


<L  +  0,  $_  )  + 

lx  lxy  ly  1 xx/ 


□fed  -  P2)[(-l6  -  10ar  +  a2)  +  (16  +  22a  +  7a2 ) p2J <^x2  - 
^(1  -  p2)a  J^(6  +  5a)  +  3(2  +  a)p2  <tL'£Ly2  - 


g(l~P  )  (2  +  a)  ^Ixx  +  ”  P2)  a(6  +  5a)  + 


(8  +  10a  +  3a2)p2  <t>lx  f[($lx2  -  $ly2)dx  +  2$lx$ly  dy 
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The  corresponding  expression  for  for  supersonic  flow  is 

o  . .  o 


obtained  by  simply  replacing  o  by  -a  with  a  = 


7  +  1 


M 


Moo2  -  1 


and 


by  _p2  with  M  ^  .  i.  The  physical  plane  variables  X  and  Y 

are  easily  inserted  into  both  equations  (l8)  and  (22)  by  means  of  the 
transformation  equations  (10), 


x  =  X 

y  =  PY 


It  is  pointed  out  that  the  forms  of  the  two  particular  integrals, 
equations  (l8)  and  (22),  derived  in  this  paper  are  identical  for  both 
subsonic  and  supersonic  flow.  The  apparent  differences  are  caused  by 
a  change  in  sign  of  the  parameter  P2.  Thus,  p2  and  a  are 

positive  for  both  subsonic  and  supersonic  flow.  Actually,  of  course, 
the  functions  represented  by  <j>^,  $£,  •  •  •  are  different  for 

the  two  types  of  flow.  For  subsonic  flow,  these  functions  are  derived 
from  analytic  and  nonanalytic  functions  of  z  and  z;  whereas  for 
supersonic  flow,  they  involve  the  real  "characteristics”  variables  x  ±  Py. 

Note  that  the  last  term  of  the  expression  on  the  right-hand  side 
of  equation  (22)  contains  the  indefinite  integral 


It  is  obvious  from  the  corresponding  complex  integral  in  equation  (21) 
that  the  integrand  of  equation  (23)  is  an  exact  differential.  This  fact 
can  also  be  easily  verified  with  the  help  of  Laplace Ts  equation, 


Thus, 


®Lxx  +  $lyy  =  0 


* 


-  v)  -  s(2*l*v) 


(24) 
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Note  that  in  general  it  is  simpler  to  perform  the  complex 

P  P 

integration  /  w^z  dz  rather  than  transform  to  a  real  integral 
then  perform  the  integration. 


and 


Langley  Aeronautical  Laboratory 
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